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16C.  LARGE-AMPLITUDE,  HIGH-RATE  ROLL  OSCILLATIONS  OF  A 65° 

DELTA  WING  AT  HIGH  INCIDENCE 

Neal  M.  Chaderjian  and  Lewis  B.  Schiff 
NASA  Ames  Research  Center 
Moffett  Field,  California,  U.S.A. 


INTRODUCTION 

The  IAR/WL  65°  delta  wing  experimental  results  provide  both  detail  pressure  measurements  and  a wide  range 
of  flow  conditions  covering  from  simple  attached  flow,  through  fully  developed  vortex  and  vortex  burst  flow,  up  to 
fully-stalled  flow  at  very  high  incidence.  Thus,  the  Computational  Unsteady  Aerodynamics  researchers  can  use  it  at 
different  level  of  validating  the  corresponding  code.  In  this  section  a range  of  CFD  results  are  provided  for  the  65° 
delta  wing  at  selected  flow  conditions.  The  time-dependent,  three-dimensional,  Reynolds-averaged,  Navier-Stokes 
(RANS)  equations  are  used  to  numerically  simulate  the  unsteady  vortical  flow.  Two  sting  angles  and  two  large- 
amplitude,  high-rate,  forced-roll  motions  and  a damped  free-to-roll  motion  are  presented.  The  free-to-roll  motion  is 
computed  by  coupling  the  time-dependent  RANS  equations  to  the  flight  dynamic  equation  of  motion.  The 
computed  results  are  compared  with  experimental  pressures,  forces,  moments  and  roll  angle  time  history.  In 
addition,  surface  and  off-surface  flow  particle  streaks  are  also  presented. 


LIST  OF  SYMBOLS  AND  DEFINITIONS 


B wing  span,  (in) 

c root  chord,  (in) 

c0  mean  aerodynamic  chord,  (in) 

Cp  pressure  coefficient  =(p-p0)/qs 

q rolling  moment  coefficient  =f/qsB 

CN  normal  force  coefficient  =N/qs 

f frequency,  (Hz) 

k reduced  frequency  =7tfB/V0 

(.  rolling  moment,  (lbs-in) 

M«  Mach  number 

m pitching  moment,  (lbs-in) 

N normal  force,  (lbs) 

n yawing  moment,  (lbs-in) 

p pressure,  (psi) 

Po  static  pressure,  (psi) 

q dynamic  pressure,  (psi) 

Re  Reynolds  number,  based  on  root  chord 

s wing  area,  (in2) 

S semi  span,  (in) 

To  static  temperature’  (°C) 

t time  (sec) 

Vo  free  stream  velocity  (ft/sec) 

x,y,z  body  axes  coordinates 

XCp  center  of  pressure  in  x axis,  (in) 

a angle  of  attack,  (°) 

CT  sting  angle  (between  body  axis  and  tunnel  axis),  (°) 

<})  roll  angle,  (°) 

<(>o  mean  roll  angle  or  initial  roll  angle,  (°) 

A<p  amplitude,  (°) 
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0 

roll  angular  rate,  (rad/sec) 

ccw 

counter  clockwise 

cw 

clockwise 

CFD 

Computational  Fluid  Dynamics 

NSS 

Navier-Stokes  Simulation 

RANS 

Reynolds-averaged,  Navier-Stokes  Equations 

FORMULARY 

General  Description  of  model 


1.1 

Designation 

1AR  Delta  Wing 

1.2 

Derivation 

IAR  Dynamic  Experimental  Model 

1.3 

Type 

Full  model 

1.4 

References 

Ref.  1 (Fig.  1) 

Model  Geometry 

2.1 

Planform 

Delta  wing-body,  See  Fig.  1 

2.2 

Aspect  ratio 

1.866 

2.3 

Mean  aerodynamic  chord 

16.323  in 

2.4 

Root  chord 

24.485  in 

2.5 

Span 

22.835  in 

2.6 

Reference  center 

13.875  aft  of  the  apex 

2.7 

Leading  edge  sweep 

65° 

2.8 

Trailing  edge  sweep 

0° 

2.9 

Taper  ratio 

0 

2.10 

Twist 

0° 

2.11 

Dihedral 

0° 

2.12 

Area  of  planform 

279.486  in2 

2.13 

Leading-edge  bevel  (leeward) 

9°  (perpendicular  to  leading-edge) 

2.14 

Leading-edge  bevel  (windward) 

9°  (perpendicular  to  leading-edge) 

2.15 

Trailing  edge  bevel  (leeward) 

9°  (perpendicular  to  trailing  edge) 

2.16 

Trailing  edge  bevel  (windward) 

9°  (perpendicular  to  trailing  edge) 

2.17 

Leading-edge  radius 

0.020  in 

2.18 

Tolerance  of  leading-edge  radius 

±10% 

2.19 

Definition  of  profiles 

0.375  inch  thick  flat-plate  wing  with  double-bevelled  (18° 
included  angle)  sharp  leading  and  trailing  edge 

2.20 

Center  body 

Ogive-cylinder 

2.21 

Form  of  wing-body  junction 

Bevelled,  see  Fig.  1 

2.22 

Form  of  wing  tip 

Sharp 

2.23 

Control  surface  details 

None 

2.24 

Center-body  diameter 

3.150  in 

2.25 

Radius  of  forebody 

r = ^24.1032  —(12.243  — x)2  —22.528  in 

2.26 

References 

Ref.  1 (Fig.  1) 
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CFD 

Grid  Details 

3.1 

RANS  grid  size 

67  axial  x 209  circumferential  x 49  normal  points  (baseline  grid); 
1 13  x 421  x 97  points  (finest  grid),  See  Fig.  2 

3.2 

Additional  Remarks 

Full-body  grids  used  in  all  cases;  zonal  grids  used  in  axial  directio 
to  fit  machine  memory;  zonal  boundaries  are  one-to-one  matching 

CFD 

Code  used 

4.1 

RANS  code 

Novier-Stokes  Simulation  (NSS)  code,  Beam-Warming,  block  or 
diagonal,  central  differencing,  blended  2nd-  and  4th-order 
dissipation,  reduced  dissipation  in  boundary  layer 

4.2 

Turbulence  model 

Baldwin-Lomax  with  Degani-Schiff  modifications,  no  fixed 
transition 

4.3 

Computational  time  step 

3.62  x 10"^  < t < 5.0  x lO'^,  6.67  x 10"®  sec  < At  < 9.0  x 10"®  sec 

4.4 

Computation  time 

80  hours  per  oscillation  cycle  on  a CRAY  C-90  single  processor  - 
block  version  (15,000  steps  per  cycle  of  oscillation) 

4.5 

Additional  remarks 

Unsteady  computation  started  with  steady  solution  at  <|>max. 
Solution  converged  after  2-3  cycles 

4.6 

Reference  on  code 

Ref.  [5] 

Model  Motion 

5.1 

Mode  of  applied  motion 

Sinusoidal  roll  oscillations  and  free-to-roll  motion  about 
longitudinal  axis  of  symmetry 

5.2 

Range  of  amplitude 

28.2°,  3 1.9°,  40.0° 

5.3 

Reduced  frequency 

f = 7 Hz,  10  Hz;  k = 0.14, 0.20 

5.4 

Additional  Remarks 

oscillations  about  <]>o  = 0.0°,  28.0° 

Boundary  Conditions 

6.1 

Mach  Number 

0.27 

6.2 

Reynolds  Number 

Rec  = 3.67  x 10®,  based  on  root  chord 

6.3 

Temperature 

300°  K 

6.4 

Range  of  model  incidence 

15°  and  30° 

6.5 

Definition  of  model  incidence 

Model  incidence  defined  relative  to  model  axis  of  symmetry 

6.6 

Additional  Remarks 

Distance  of  far  field  boundary  is  5 root  chords  normal  to  wing,  2 
root  chords  upstream  and  downstream  of  wing 

Data 

Presentation 

7.1 

Static  Cases 

Roll  moment  versus  roll  angle 
Normal  force  versus  roll  angle 
Center  of  pressure  versus  roll  angle 

Leeward  surface  pressure  distributions  for  various  roll  angles 
Surface  flow  patterns  for  various  roll  angles 
Vortex  breakdown  point  versus  roll  angle 

7.2 

Forced  Roll  Oscillations 

Instantaneous  roll  moment  versus  roll  angle 
Instantaneous  normal  force  versus  roll  angle 
Center  of  pressure  versus  roll  angle 

7.3 

Free-to-Roll  Oscillations 

Roll  angle  history  versus  time 

7.4 

Sample  illustrations 

Fig.  3 to  Fig.  10 
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7.5  Additional  Remarks 


All  above  quantities  are  compared  against  IAR  experiments  (Case 
No.  2,  No.  3,  No.  5,  No.  6 and  No.  10) 


Personal  contact  for  further  information 

Dr.  Neal  M.  Chadeijian,  T27B-2,  NASA  Ames  Research  Center,  Moffett  Field,  California  U.  S.  A 
Phone:  (650)  604-4472,  E-mail:  nchaderjian@mail.arc.nasa.gov 
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Fig.  3 Effects  of  grid  refinement  and  zonal  boundary 
condition  treatment  on  the  pressure  coefficients 
M„=0.27,  a=15°,  6=0,  Re=3.67  million 


Fig.  4 Comparison  of  computational  and  experimental  rolling 
moment  coefficients  for  dynamic  and  static  cases 
M„=0.27, 0=15°,  A<t>=40°,k=0.14,  Re=3.67  million 
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Fig.  5 Comparison  of  mean  computed  and  experimental 
rolling  moment  coefficients  for  static  roll  angles 
M„=0.27,  O=30°,  Re=3.67  million 


Fig.  6 Dynamic  and  static  rolling-moment  coefficients 
M„=0.27, 0=30°,  $,,=0”.  A(|*=28.2°,k=0.20,  Re=3.67  million 


Fig.  7 Dynamic  and  static  rolling-moment  coefficients 
M„=0.27,  0=30”,  (Jj0=28°,  A<|>=31.9°,  Re=3.67  million 


Fig.  8 Time  history  of  roll  angle  for  free-to-roll  motion 
M„=0.27, 0=30°,  <t)o=40.5°,  Re=3.67  million 


